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Abstract. The self-consistent harmonic approximation is extended in order to 
account for the existence of Klein factors in bosonized Hamiltonians. This is 
important for the study of finite systems where Klein factors cannot be ignored 
a priori. As a test we apply the method to interacting spinless fermions with 
modulated hopping. We calculate the finite-size corrections to the energy gap 
and the Drude weight and compare our results with the exact solution for special 
values of the model parameters. 



PACS numbers: 71.10.Pm, 71.30.+h 

1. Introduction 

Bosonization techniques have been extensively used in the study of one-dimensional 
electron and spin systems |2 El • The success of the method is based on the fact 
that the low energy properties of fermions are determined by the states close to the 
Fermi surface which in one dimension consists only of two points. When the spectrum 
is linearized around the Fermi points the Hamiltonian can be expressed in terms 
of bosonic operators associated with particle-hole excitations of momentum q. In 
addition, scattering processes between left and right Fermi points require to introduce 
extra operators, the so-called Klein factors. 

Perturbations like impurity scattering or a modulation of the hopping lead to 
nonlinearities in the bosonized Hamiltonian, for which an exact solution is known 
only in some special cases In general one has to resort to approximative methods 
like renormalization group calculations. Another more intuitive method is the self- 
consistent harmonic approximation (SCHA) where the non-linear terms are replaced 
by a harmonic potential with parameters to be determined self-consistently according 
to a variational principle for the energy. 

The SCHA has been successfully applied to various non-linear field theories 
|Kll?)l[7llHII??l 1101 ITT] . In the context of bosonized Hamiltonians, however, the existence 
of Klein factors has been ignored in these approaches. This may be justified for infinite 
systems (3|, but in general the Klein factors have to be treated carefully as pointed out 
for example in the context of impurity models and two-leg ladders ^1 ^] . In 
the following we present an extension of the SCHA which treats the bosonic fields and 
the Klein factors on equal footing. As a prototype model for this study we consider the 
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lattice model of spinless fermions at half filling including nearest-neighbor interaction 
and a modulation of the hopping due to a periodic lattice distortion. 

In the next section the Hamiltonian of this model is bosonized and a trial 
Hamiltonian is constructed which serves as the basis for the SCHA. In Sect. 3 we 
calculate the energy gap and the Drude weight and compare our results with exact 
results that are available for certain values of the model parameters. Finally we 
conclude with a brief summary in Sect. 4. 



2. Model and formalism 



We consider a one-dimensional model of spinless fermions with static dimcrization 

i 

-M/^niTii+i (1) 

i 

where u is the dimerization parameter that leads to a periodic modulation of the 
hopping amplitude t, and V is the strength of the nearest-neighbor interaction. At 
half filling and in the range —2 < V/t < 2 the bosonized Hamiltonian of this model 
reads |5] 



FiFn e^'^(^) + h.c. (2) 







where (j) and 9 are conjugate fields. The renormalized Fermi velocity v = 
TTt sin(2ry) / (77—277) and the Luttinger parameter g = n/Ari are related to the interaction 
according to = — 2^003(277). The normal ordering operation (: :) is introduced 
to avoid the divergences associated with unphysical states assumed by the Luttinger 
model. The parameter a is needed in the bosonization formalism to remove ultra-violet 
divergences in certain fc-summations, and appears in Eq. |(2J) since the exponential 
exp{2i(f>{x)) is not normal ordered. For a lattice model, 1/a can be interpreted as 
a kind of effective band-width 2 , whereas in the continuum the limit a — > has 
to be performed in order to obtain the correct 5-function anticommutators of the 
fermionic field operators pp. N and J are the charge and current operators, defined as 
= Nl -\- Nn and J = — N^, respectively, where Nl and Njf count the number of 
left and right moving particles with respect to the filled Fermi sea. The Klein factors 
^L/R ^^'^ ^L/R I'ise or lower the number of left and right moving fermions by one 
(which no combination of bosonic operators can ever do) , and they assure also that 
fermions of the two different species anticommute. The dimerization, described by the 
non-linear term of the Hamiltonian, breaks the conservation of the current J and is 
responsible for the opening of an energy gap. 

To discuss the gap formation quantitatively, we use the self-consistent harmonic 
approximation (SCHA) which has originally been introduced for the sine-Gordon 
model [H]. The idea is to construct an exactly solvable trial Hamiltonian iJtr where 
the bosonic fields in Q are decoupled from the Klein factors and where the non- 
linear terms ~ e^^*"^ are replaced with a quadratic form ~ (jp' . Accordingly we chose 
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Ht, = + Hf^ as the sum of two commuting parts 

i^tt - rp {-{d.<l>f+v9{d.er + — j (3) 
Jo 27r [5 vg J 

and 

= -^B{FtFn - F^F,) + ^ (4) 

each of them depending on a single variational parameter, A and B, respectively. 
Notice that from the standard trial Hamiltonian alone (Eq. only the trivial 

solution A = can be obtained. Therefore the non-standard term Eq. is necessarily 
required to describe the phase transition in the system. Since we are only interested 
in properties at zero temperature we use the normalized ground state |V'o) of iJtr as 
a variational wave-function for the Hamiltonian (0). This provides us with an upper 
bound E for the ground state energy E oi ^ due to the inequality 

E<E^{^o\H\i>o)- (5) 

Minimizing E with respect to the variational parameters A and B yields the following 
set of equations: 

A2 = -(4i;.giu/a)£;^(B)e-2(<#''>tr 

(6) 

B ^ (tu/7ra)Le-2(-^ 

where Eq{B) is the ground state energy of and Eq{B) — dEQ{B) /dB. Since 
is bilinear in the field operators it is straightforward to calculate the equal time 
correlation function of the phase field (f>{x) entering Eq. 

where the /c- values in the sum arc multiples of 2tt/L. 
3. Energy gap and Drude weight 

In order to solve the set of equations we have to calculate the ground state energy 
Eo{B) of the trial Hamiltonian H^. The products of Klein factors F^Ffi and F^F^ 
appearing in change the quantum number of the current operator from J to J±2. 
iJj^ can therefore be regarded as a particle moving in a harmonic potential along the 
J-axis. For L ^ 00 the potential energy term in vanishes and the kinetic energy 
term with imaginary hopping iB yields the ground state energy Eq{B) = —2B. In this 
limit we reproduce the results of Ref. ^I], where the Klein factors have been ignored. 
Replacing the sum in Eq. (jT)) by an integral we obtain 

('^^-fln^ (8) 

where Aq — 2ve~^ /a, 7 = 0.5772 is Eulcr's constant, and A ^ Aq has been assumed. 
Inserting Eq. ||SJ| into Eq. © and using E'q{B) = —2 yields 

Ao V"o/ 

with uq = e~'' Ao/Agt. For g > 2 Eq. diverges as u ^ and there is only the 
trivial solution A = 0, i.e. the line g = 2 marks the transition from a gapless to a 
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gapped phase in the g-u plane. In the spinless ferniion model this corresponds to 
the line V/t = —a/2 in the V-u phase diagram. The value gc = 2 obtained within 
the variational approach is in accordance with the exact result based on a mapping 
of the bosonized Hamiltonian on the Ashkin- Teller model for w — > 0. However, 
according to renormalization group calculations |16| and confirmed numerically using 
density matrix renormalization group '11' , the phase boundary in the g-u plane shifts 
to larger values of g with increasing u while the SCHA gives a vertical line. On the 
other hand, the exponent l/{2—g) characterizing the opening of the gap is exact [S]. In 
particular for non-interacting fermions (V — 0) the Luttinger liquid parameter is g = 1 
and therefore A oc m in agreement with the exact energy gap Aox = 4tM. At V = 2t, 
where Umklapp scattering becomes relevant and leads to the opening of a correlation 
gap, Eq. (jnj yields A esc u^^^ which agrees with the exact result Aox oc v?!'^ j \J\ lnu| 
up to the logarithmic correction |17j . 

In the following we consider a system of finite length L. In order to calculate the 
ground state energy of it is convenient to switch to the momentum representation, 
J — id/dpj, where the trial Hamiltonian reads 

with periodic boundary conditions for the wave-function, ^'(pj -I- tt) = ^'(pj). The 
corresponding Schrodinger equation is known as Mathieu's equation. Accordingly, the 
ground state energy of Hf^ has the following asymptotical behavior j2J 

E,{B)^{ f ^^«^ (11) 

' \ -2B iovLB:^v ^ ' 

Combining the two equations © one obtains 

A^L^ = -4:nvgE'Q{B)BL (12) 

which can be used to relate the two limiting cases of Eq. Hll|) with two different 
physical situations: The condition LB ^ v is equivalent with LA ^ v and therefore 
the results for ((/)^)tr and A(w) are the same as for the infinite system i.e. the energy 
gap depends algebraically on the dimerization. On the other hand, the condition 
LB <C V implies that LA ^ v. In this case the energy gap is purely due to the finite 
system size and of order 2ttv/L. Thus the crossover from a finite-size gap to a true 
dimerization gap coincides with the crossover between the regions where Klein factors 
are relevant or can be ignored. 

In general, the set of equations © can only be solved numerically. Fig. 1 shows the 
size-dependent energy gap A(L) = y/ {2nv/L)'^ -\- A^, i.e. the lowest excitation energy 
of the bosonic modes, as function of u for several values of the system size L; notice 
that we consider a system with periodic boundary conditions, where the smallest wave 
number is 2n/L. The Luttinger parameter is g = 3/4 which corresponds to V/t — 1. 

We now turn to the calculation of the Drude weight D within the SCHA 
formalism. The Drude weight or charge stiffness is the central quantity to characterize 
charge transport at zero frequency. At T = the real part of the electrical conductivity 
at frequency oj is of the form — 2ttDS{uj) -\-aj:cg{Lu) with crrog(w) — > for lo —^ in 
a system without impurities. Therefore D — characterizes an insulator while D > 
describes an (ideal) conductor. The simplest way to calculate D is due to Kohn [TH) 
who has shown that the Drude weight can be expressed as 



2 d^2 



(13) 

ip=0 
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Figure 1. Energy gap A(L) (in units of t) as function of the dimerization 
parameter m for L = 100, 200, 500, 1000, 2000 (from top to bottom). The dashed 
hne is the analytic result of Eq. j^J valid for L —> oo and A "C Aq . 



where E(ip) is the ground state energy of a ring of circumference L which is threaded 
by the flux (/?. Alternatively, the parameter (p can also be associated with a change 
of boundary conditions, i.e. ip = corresponds to periodic and ip = n to antiperiodic 
boundary conditions. In the fermionic model the hopping parameter t picks up 
a phase factor exp(±i(p) when a flux ip is applied. In the bosonized form (O the 
only modification is that the current operator J has to be replaced by J + 7. 
Correspondingly, we modify the i?-dependent part of the the trial Hamiltonian and 
write 

, 2 



TTvg 



(14) 



Applying the same procedure as in the case p = yields a variational estimate E for 
the ground state energy which now depends on the flux p, i.e. E = £'(A, B, p) where 
A and B, obtained from the gap equations 0, are also functions of p. On first sight 
it might seem impossible to calculate the Drude weight using Eq. (|13|) since there is 
no analytical solution of the gap equations © even for 93 = 0. However, a closer look 
reveals a great deal of simplification. From the second derivative 



d^E d^E 



dpi'^ 



dA 
dE d^B 



d^E 



dB^ 



d^E 



dBdp^ dp^ 



one retains only 



d^E 



dp^ 



d^E{A,B,p) 



ip=a 



dB 
dp 



^ dE d^A 

+ 



OA dp^ 



mixed terms 



(15) 



(16) 



since i) ^ — ^ ~ due to the minimum condition of the energy and ii) 
at (/? = due to symmetry (A and B are even functions of p). Since 
the trial Hamiltonian Htr = + H^{p) consists of two commuting parts we may 



dA 

dip 



OB 

dp 
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write £^(A, B, (p) = £'(A) + Eq{B, ip) where £;(A) depends only on A and Eo{B, ip) is 
the ground state energy of H^.{ip) i.e. we obtain the simple result 



2 dip 



2 



(17) 

where B is given by Eq. I^. To proceed we again represent H^.{p) in terms of the 
Mathieu equation H10|l where now the boundary conditions are ^'(pj + 7r) = e*'^\I'(pj). 
It is now straightforward to calculate the Drude weight for a finite system of size L in 
the gapped phase g < 2. In the "finite size gap" region (LA ^ v) we obtain 

D = Doil-^ + ...) (18) 

where q = ^LB/irvg oc uL^"^ and Dq — vg/2n is the Drude weight of an unperturbed 
Luttinger liquid. In the opposite limit (LA ^ v) which corresponds to a Mathieu 
equation with a large cosine potential the variation of the ground state energy with 
change of boundary conditions is exponentially small (as expected from the WKB 
approximation). A more careful treatment j21j yields 

D^D^-(—X'\-^^^'-^<^ (19) 
t: \vg J 

where A and u are related via Eq. I^. In the special case of free fermions i.e. for 
(7 = 1 and v = 2t we may compare the results obtained within the SCHA with the 
exact Drude weight (see Appendix) 

D,.c.i?o(^^)'^^-^'^-/^^ (20) 

Although the expected exponential behavior is recovered in the SCHA result, the 
different numerical factor in the exponent shows the limitations of the method. 



4. Summary 

We have constructed an extension of the SCHA in order to account for the existence 
of Klein factors in bosonized Hamiltonians with non-linear perturbations. As an 
application we have investigated a model of spinlcss fermions with modulated hopping. 
For the infinite system, both the value of the Luttinger parameter gc — 2 where the 
transition from a gapless to a gapped phase takes places and the exponent 1/(2—17) 
that characterizes the opening of the gap for u — s- are correctly obtained within 
the SCHA. However, the bending of the phase boundary for finite values of u is not 
reproduced. When considering a finite system Klein factors cannot be ignored a priori. 
Within our approach it turns out, that the crossover region from a finite size gap to 
a true dimerization gap coincides with the crossover to the region where the Klein 
factors become relevant. 

The Drude weight reflects the sensitivity of the system with respect to a change 
of boundary conditions and is related to the properties of the current operator J in 
the bosonized version of the Hamiltonian. In a finite system with an energy gap A the 
Drude weight is expected to be nonzero but exponentially small, D ~ exp(— const-LA). 
Our extended version of the SCHA allows to calculate the Drude weight in the 
insulating phase and we confirm the exponential behavior. 

Using the same concepts the method can also be applied to study finite size effects 
in more complex models with nontrivial phase diagrams. 
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Appendix A. 

For V = the ground state energy of Hamiltonian (1) for a system of TV = AM + 2 
lattice sites is given by 

M 



n=-M 



£;o((p) = -2t y d cos^ kn + V? sur kn (A.l) 



where fc„ = {2'im-\-if)/L and the lattice constant is set to one. We introduce z„ = e**''" 
and the function /(z) = \/{z^ — e"^) which has single poles at z„ with residua 
Res/(2;)|z=z„ = Zne~^'^ /L. Expressing Eq. (|A.ip in terms of z„ we may replace 
the sum by a contour integral and obtain 



^0 <^ = ^ <P TT- , T A.2 

where 7 = (1 + — u^) and C is a contour that encloses the singularities of /(z). 

We chose C to be composed of two circles with radii slightly larger or smaller than one, 
respectively. Substituting z — > 1/z the integral along the inner circle can be mapped 
on the integral along the outer circle. The square root in the numerator has branch 
cuts along the imaginary axis in the interval [— 2/i,yi] and for \y\ > yo where ±iyo 
and ±iyi are the zeroes of the function under the square root. Now we deform the 
integration contour along the branch cut from iyo to ioo and obtain for the energy 
difference AE = Eoin) - Ea{0) 



Jyo y y^-y-^ 



Substituting y — yields 



AE ^ l^^^nZ r dx V^^o^^i^-}-^^ (A.4) 



with xq — Arcosh7 = ln((l + u)/{l — u)). For Lxo ^ 1 the integral is rapidly cut off 
by the exponential and we may expand the square root around x = xq. The Drude 
weight is then 



^ LAE ^ 2uL ( L\+u\ ,^ , 

i?=^c.y_exp(^--ln^j (A.5) 

Expanding the logarithm for m ^ 1 and inserting Acx — 4:tu yields Eq. H20|) . 
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